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Energy level statistics of Hermitian random matrices H with Gaussian indepen- 
dent random entries IIi>j is studied for a generic ensemble of almost diagonal random 
matrices with ~ 1 and = bj^{\i — j\) <^ 1. We perform a reg- 

ular expansion of the spectral form-factor K{t) = 1 -|- 6i^i(r) -|- b'^ K2{t) + . . . in 
powers of 5 <C 1 with the coefficients Km{T) that take into account interaction of 
{m + 1) energy levels. To calculate KmiT), we develop a diagrammatic technique 
which is based on the Trotter formula and on the combinatorial problem of graph 
edges coloring with {m + 1) colors. Expressions for Ki{t) and K2{t) in terms of 
infinite series are found for a generic function J^{\i — j\) in the Gaussian Orthogo- 
nal Ensemble (GOE), the Gaussian Unitary Ensemble (GUE) and in the crossover 
between them (the almost unitary Gaussian ensemble). The Rosenzweig- Porter and 
power-law banded matrix ensembles are considered as examples. 

PACS numbers: 71.23.-k, 71.23.Aii, 71.30.-hh, 02.10.Yn 



I. INTRODUCTION 



Random Matrix Theory (RMT) has proven to be a universal formalism to describe a 
great variety of complex systems ranging from nuclei to mesoscopic quantum dots^ to chaotic 
systems^. Of particularly wide application is the Wigner-Dyson RMT^^i^iSi. It is the statistical 
theory of eigenvalues and eigenf unctions of a random Hermitian matrix H whose entries 
Hi^j fluctuate as independent Gaussian random variables with zero mean (Hij) = and a 



constant variance {\Hij\'^) = const. There are three Dyson symmetry classes labeled by /3: 
orthogonal {(3 = 1), unitary (/5 = 2) and symplectic {(3 = 4) that correspond to real, generic 
complex, and real quaternionic Hermitian matrices. This theory is extremely successful in 
describing the spectrum of complex nuclei and the statistics of various observable quantities 
in mesoscopic quantum dots^^. 

The key property of the Wigner-Dyson RMT is that the variance does not change 

with increasing the distance \i — j\ from the diagonal. In order to illustrate the properties 
of eigenvectors one can invoke one-dimensional chain of sites 1 < i < (A^ is the matrix 
size) with on-site energies Ha and hopping matrix elements Hij. Then the above property 
of off-diagonal entries implies the possibility of hopping throughout the entire chain which 
results in the delocalized character of eigenvectors. If on the contrary the off-diagonal matrix 
elements Hij are nonzero only inside the band |« — j| < A, all the eigenvectors turn out to 
be localized in the same way as in a quasi one-dimensional disordered wire^. This banded 
random matrix (BRM) theory is also relevant for chaotic systems under a time-periodic 
perturbation^. 

One well-known extension of the Wigner-Dyson RMT is the model of Rosenzweig and 
Porter^ (RPRM) with the variance = B'^/N^"' that does not depend on the distance 

but depends on the matrix size. The spectral statistics of RPRM at different parameter 
a can range from almost uncorrelated levels (if a > 1 ) to a strong level repulsion which 
is similar to a weakly perturbed Wigner-Dyson ensemble (if a = 1/2), see, for instance, 
the papers [jlfl^il^l l and references therein. Exact calculation of the two-levels correlation 
function of RPRM with the complex hopping entries Hi^j and a = 1 shows that the 
spectral statistics falls neither into the Poissonian nor into the Wigner-Dyson universality 
classesiS. It is often called "the regime of crossover" . Much less is known about RPRM with 
either real or quaternionic hopping elementsi^. 

Without lost of generality, RPRM can be presented by an ensemble of the matrices 
Hrp = a + {B/N°')B with a superposition of a diagonal matrix A and a full Wigner- 
Dyson matrix B of any symmetry (/? = 1,2,4). If A is the random matrix possessing 
the Poissonian level statistics and E = Q the model describes spectral properties of a 
classical integrable system. By increasing B one can explore a transition from integrable to 
a classically chaotic system with or without the time-reversal symmetrj*i^. 

Recently another Gaussian RMT that interpolates between the Wigner-Dyson RMT and 



BRM theory attracted a considerable attention. This is the power law banded random matrix 
(PLBRM) theorj*i^*i^»i^*ii for which is nearly constant inside the band < A and 

decreases as a power-law function ~ l/|z— for > A. The special case a = 1 

is relevant for description of critical systems with multifractal eigenstate o^'^'^^'^i'^^i'^'i'i'^^i'^^ , in 
particular for systems at the Anderson localization-delocalization transition point. The case 
a > 1 corresponds to the power-law localization which can be found in certain periodically 
driven quantum-mechanical systems^S. 

The progress in BRM and PLBRM theories became possible because of mapping2*i^ onto 
the nonlinear supersymmetric sigma-model^^ that allowed to obtain rigorous results by using 
various powerful methods of the field theory. However, such mapping is only justified if the 
bandwidth A ^ 1. In the opposite case where all the off-diagonal matrix elements are 
parametrically small compared to the diagonal ones, no field-theoretical approach is known 
so far. Yet such almost diagonal RMT may possess nontrivial properties because of the slow 
decay of the off-diagonal matrix elements {\Hij\'^) with increasing |« — j|. For instance it is of 
fundamental interest to study the spectral statistics in systems with power-law localization 
that takes place in the power-law banded random matrix ensembles at a > 1. Another 
problem to study is the critical almost diagonal PLBRM. It is known that the eigenvectors 
of PLBRM with a = 1 remain multifractal for an arbitrary small value of A |l9|. This 
means that the typical eigenfunction is extended though very sparse at small A. Thus 
almost diagonal random matrices may display the localization-delocalization transition with 
changing the exponent a as well as their large bandwidth counterpart^^. This transition 
has been studied numerically-- for A ~ 1 but little is known about it in the limit of almost 
diagonal random matrices. 

The goal of this paper is to develop a formalism that would allow to describe the spectral 
statistics for a wide class of almost diagonal random matrices with a generic behaviour of 
= 6^jF(|i — j|) and ~ 1. The parameter b may or may not depend on the 

matrix size A^. The principle requirement is that 6 ^ 1 is small. It controls the smallness of 
the off-diagonal matrix elements and will be used as an expansion parameter for the spectral 
correlation functions. 

A natural way to proceed with such an expansion is to develop a perturbation theory in 
the parameter b starting from the diagonal matrix as a zero-order approximation. However, 
a naive expansion of the spectral form-factor K{t) up to a finite order in b may fail, since 



the parameter b enters in the product bNr with the matrix size N and the time r 
measured in units of the Heisenberg time. In this combination, the small parameter b 
can be compensated by the (large) product Nr. In order to be able to obtain the regular 
expansion in powers of b, K{t) = 1 + bKi{T) + b'^K2{T) + . . ., in the thermodynamic limit 

^ cxD at arbitrary Nt, we develop a diagrammatic technique that gives Km{T) in a form 
of infinite series in powers of bNr. It is essentially a kind of locator expansion^^ adjusted 
to the problem of spectral statistics. It formalizes and generalizes the idea of Ref. 1^ where 
the resonance pairs of states i,j with iHa — Hjj\ < \Hij\ have been considered to find the 
statistics of multifractal eigenstates. In our approach, the coefficients Km{T) result from an 
interaction (via the off-diagonal elements Hij) of m + 1 energy levels. We call this expansion 
a virial expansion by analogy with the expansion of thermodynamic functions of a dilute 
system in powers of density with virial coefficients that take into account collisions of m + 1 
particles. We stress that what we are doing below is not an approximation of two, three, or 
more resonant levels. It is rather the classification of exact perturbation series by a number 
of the interacting levels involved. The question of a contribution from non-resonant levels 
does not arise, since at a given number of the interacting levels all possible relations between 
\Hii — Hjj\ and \Hij\ are taken into account. 

The most serious problems on that route are i) dealing with non-commutative matrices 
= diag{Hii] and V = H — H^; ii) combinatorial coefficients in the infinite series 
in bNr. The first problem is solved by using the Trotter formula (see Section ITTT|l . The 
second problem reduces to a particular case of the general problem of graph coloring that is 
amenable to an exact solution (see Section HVl and IVjl . 

As a result we present the series for /^^(r) with combinatorial coefficients that do not 
depend on the particular choice of jF(|z — j\) in the variance = 6^jF(|z — They 

can thus be applied to any Gaussian ensemble of almost diagonal random matrices. The 
particularly simple expressions are found for m = 1, 2 which correspond to two and three 
interacting levels. 

The general formulae obtained in this way are illustrated in Section IVII Al for RPRM 
with b = B/N; T{\i - j\) = 1 and in Section IvITbI for a critical PLBRM with b = 
B; J-'{\i — j\) ~ l/\i — Jp- In both examples, the constant B will be taken small: i3 ^ 1. 
As a result, the correction to the Poisson behaviour K{t) = 1 has been obtained in the 
limit of infinite matrix size N ^ oo. The existence of such a correction can be interpreted 



as a rigorous detection of delocalization in the framework of a formalism that starts from the 
diagonal random matrix where all states are localized. Thus our theory offers a controllable 
way of obtaining delocalization from localization. This is just complementary to conventional 
approaches based on the supersymmetric sigma-model^i that always starts from delocalized 
(e.g.diffusive) modes. 

II. BASIC DEFINITIONS 

Let us consider a Hermitian RM of size N x N , N ^ 1, from a Gaussian ensemble. We 
assume that entries of the matrix are random and independent. The RM is the Hamiltonian 
H of the matrix Schrodinger equation 

where e„ and ipn are the eigenvalues and eigenvectors, respectively. Define statistical prop- 
erties of the matrix entries: 

{H,,) = 0- {Hl) = -^- {\H,,\^)=b''n\^-J\).^^J■. (1) 

where J^{\i — j\) > Q is a smooth function of its argument, and the parameter b is small: 

6< 1. 

The condition 6 <^ 1 means that RM is almost diagonal. The parameter 13 corresponds 
to the Dyson symmetry classes: (3goe = 1, Pgue = 2. The brackets (...) denote the 
ensemble averaging. If the ensemble is Gaussian the averaging of a function F{A) over a 
random variable A reads: 

We concentrate on the level statistics which is characterized by the density of states 
p{E) = J2n ^{E — e„) and its multi-point correlation functions. For example, the two-level 
correlation function R{uj) is defined as: 

fl(,), ((p(W^2y2))) ^ {{ai»-(ai>-(a)(i>. (2) 



The Fourier transform of R{u) is known as the spectral form-factor K{t): 

K{t) = / e''"*i?(cj)dcj. (3) 



We will see below that for almost diagonal RM the representation of spectral statistics in 
the time domain turns out to be convenient. Therefore from now on we use the form factor 
instead of the correlation function. We rescale time by the mean level spacing A = 1/ ( p(0) ) 
introducing the dimensionless time t = t A. Our goal is to develop a regular perturbative 
expansion for K{t) in powers of small parameter 6^1. In the limit of small time the 
spectral form-factor K{t 0) is linked to the other important spectral characteristics 
called the level compressihiliti^: 

X = lim( lim K{t)) . (4) 

Let us take a window of the width bE^ 5E//S. = n <^ N , m the energy space centered 
at = and calculate the number n of levels inside the window at some realization of 
disorder. The level number variance is Tj2{n) = {{n — n)^). The level compressibility is by 
definition the limit 

X = Jim ( lim ^] . (5) 

The level compressibility contains an information about the localization transition: x ranges 
from XWD = for the Wigner-Dyson statistics with extended wave functions and a strong 
levels repulsion to xp = 1 in the case of localized wave functions and uncorrelated levels 
with the Poissonian distribution. The intermediate situation with < Xcrit < 1 is inherent 
for the critical regime of multifractal wave functions^i. 

The main object of our further analysis is the following correlation function in the time 
domain: 

KN{r) = ^(( Tr e-'^^/^ Tr e'^^/^ )) = i^o(r) + bK.ir) + h^K^ir) + ... (6) 

For the constant mean density of states K{t) coincides with K{t). However, they are 
different if {p{E)) essentially depends on energy E. In analogy with Eq.(jH) one can define 
the quantity xo = liniT-»o( liniAf-^oo -^Af(T)) . It turns out that at small b there is a simple 
approximate relationship between x and Xo (see Appendix : 

1 - Xo ^ A 



V ~ 



/\ r-\-oo 



The mean density of states for the Gaussian ensemble of almost diagonal RMs with either 
localized or (sparse) fractal eigenstates is close to the Gaussian distribution of the diagonal 

Thus with an accuracy of 0{b) the unfolding factor in Eq.(|ZI) can be taken as ~ a/2. 



III. THE METHOD 
A. The Trotter formula 



As far as we investigate the properties of almost diagonal RMs the Hamiltonian can be 
naturally divided into a diagonal part and a matrix of hopping elements V: 

H = H, + V. (9) 

It follows from the definition (0) that the hopping elements Hij = Vij ~ b are small 
compared to the diagonal ones iJj.j = ~ 1. However, matrices He and V do not commute 
with each other and a direct expansion of the exponentials e^'*^^'""'"^-'^ in Eq.© in terms of 
V involves serious difficulties. One possible way to overcome these problems is to represent 
g±i(H£+y)^ as a product of exponential functions containing matrices and V separately. 
To do this we use a generalization of the identity exp(a + b) = exp(a) exp(6) for non- 
commuting variables known as the the Trotter formulc^: 

n 

e^+^ = hm fe^/'^e^/") e^^^i = lim TT (e^'^^±^e^'^^) . (10) 

ji— >oo V / ?i— >oo V / 

p=l 

The Trotter formula is exact for the finite matrices A and B. Therefore, correct order of 
limits reads: 

1. 72 oo ; 

< 2. oo; (11) 

3. (if necessary) r ^ . 
In what follows we will always imply precisely this order of doing limits. 

Eq. p()|l allows to make a regular expansion in powers V . The price for that is the infinite 
product in Eq. ffTUI) . We need a proper selection rule to extract from that product terms of 
the order of 0(6°), 0(6^), Oib^), . . .. 



B. Lowest-order terms 



For a strictly diagonal matrix the spectral form-factor K{t) is calculated straightfor- 
wardly: 

^(r)lv=o = ^((Tre-^^-/^ Tre^^-/^)) = ^ E [l " (e"^^-/^)^] = 1 - e-<^^>(^/^)^ (12) 

For all Gaussian RM of the type Eq.(p5 the inverse mean level spacing increases with the 
matrix size N. The inverse mean level spacing 1/A = N is proportional to the matrix size: 
N N. For finite the normalization sum rule24 requires K{0) = 0. If however the 
limit — > oo is done prior to the limit r ^ (see Eq. ljlH) ) the normalization sum rule is 
violated and we recover the Poisson statistics result K{t) = 1. 

In order to obtain corrections to K{t) proportional to 6^ one has to expand 6=*=*^^^/" in 
powers of V in the infinite product in the r.h.s. of Eq. ljlOj) : 

f Nt \ f Nt A f Nt \ f Nt A 

. . . X exp [ ±i— 1 exp ( ±i— ^ 1 x exp ( ±i—He 1 exp ( ±i— ^ 1 x . . . . (13) 



n pairs of exponentials 

and then to perform the Gaussian averaging over and V. Note that the Gaussian average 
is zero for all terms in r.h.s of Eq.© which contain an odd number of off-diagonal matrices 
V. That is why only terms with even powers of b may result from such an expansion. For 
the purpose of explaining the details of the new formalism based on the Trotter formula we 
show how the term oc 6^ arises from it. 



1. Gaussian averaging 

Each Tre^^-^^^ in Eq.© can be represented using the Trotter formula as Tr of the 
infinite product Eq. pHj) . Therefore we have to expand either two exponentials e^*^^"^/" in 
Eq. (fT!?|l up to V or one single exponential up to V"^ setting = in all the other exponentials 
^±tVNT/n^ We will show below that only the first option survives the limit n ^ oo. 

Expanding two of exp {±i^v) in Eq.jni) up to we obtainSi: 



exp ±1 i/g ±1 V exp ±1 i/g ±1 V exp ±i ~ o 

V n \ n I \ n \ n \ n / 



(14) 



FIG. 1: Graphic illustration of Eg . (fT6|l before (a) and after (b) averaging. 



where pi,2,3 are the numbers of successive exponentials exp [±i^He) 'fused' together after 
we set V = in all exp ^±i^V^^ but two. We call ps the Trotter numbers. They obey the 
obvious restriction 

Pl+P2+P3=n, Pi,2 > 1 , P3 > . (15) 

Note that both off-diagonal matrices V in Eq. ()14|) must belong to the same trace. If 
they belong to the different traces the result is zero because Tr {exp (±i^4) V} = 0. 

In order to find the term C oc 0(6^) in Eq.® we substitute Eq. (fTl|l into Eq.ijHJ and 
perform the Gaussian averaging in accordance with Eq.((H): 



2 

C = -^(^) E(( E e-^"^-V,,e-'^^^V^,,e-^"^-xe'^-'"))U=.or„.=,(16) 

V ^ / {ps} i,j,m=l 



4 /Nt\ " " 



{Ps} 

n n 

where E means summation over the Trotter numbers E ^n-(pi+P2+P3) ^^^1 we denote 

{ps} Pl,P2,P3 = l 



N 

I / ~ \ / _ I 

I^jv(2/) = exp 



2(3 



1 ^ k 
'^N{k)^-Y.{H\^-3\)) ■ (17) 



The function TZ^lk) depends on the particular form of decay of the off-diagonal matrix 
elements in the Gaussian ensemble (see Eq.(^) and contains the sum over matrix indices 
which will be referred to below as the summation over the 'real space'. 

It is instructive to give a graphic representation of Eq. p6|) before and after averaging. In 
FiglUeach circle denotes the trace. The vertices Vij and Vj^i in the left circle are uncoupled 



before averaging (Fig^). The shadowed box indicates starting point of the Trotter chain 
Eq. lfT^ . The circle is thus divided into segments (below referred to as the Trotter segments) 

. Ntpi . Ntp2 • Ntp^ 

corresponding to e~^~^'^^ , e~^~^'^^ , and e~^~^^\ We mark the length of each segment 
by the Trotter number ps. The right circle stands for e'^^^™. The averaging over V connects 
uncoupled ^-vertices by the correlation function \Vij\'^ = 6^jF(|z — j\) denoted by a double 
line while the averaging over couples Em to either Si or Ej (FigOD). 

2. Integral over Trotter variables 

One has to sum over the Trotter numbers in Eq. p6|) prior to doing the limit n — > oo. 
To this end we introduce the Trotter variables Us = Ps/n which in the limit n — ^ oo can be 
considered as continuous. Then the summation over the Trotter numbers Ps can be replaced 
by integration over the Trotter variables. Switching from the Kronecker-delta 5p to the 
Dirac (5-function n5p 6{y) we arrive at: 

C = -4(6iVr)^ 7^iv(l)X; 

^ = JJJq d?/id?/2d?/3 5(l -yi-y2- 2/3){'Piv (1 -yi - 2/3) "^iv (2/2) + (18) 

+Vn (2/1 + 2/3)^^^(1-2/2)}. 

It is remarkable that as a result of this transformation, which is exact in n ^ 00 limit, all 
the n-dependent factors are absorbed by 2/s- This is because the number of independent 
Trotter numbers Ps in Eq. fllfjj) coincides with the number of matrices V, each coming with 
the factor 1/n. Note that this is not the case if only one exponential e±*^^^/" in Eg. ffT^ 
is expanded up to V"^ order. Then the number of Trotter summations is less by one which 
leaves an uncompensated factor 1/n after switching to the Trotter variables. Thus we 
arrive at an important conclusion that in all orders of V^-expansion one should retain only 
linear terms of expansion of the exponential e^*^^^/" 1 ± iVNr/n 

The integral Eq. lfT^ . which will be referred to as 'the Trotter integral', can be simplified 
by observing that the (5-function imposes the same arguments in both Vn functions in the 
corresponding products. Thus functions can be fused Vn{x) Vn{x) = [V2 x] and 
the Trotter integral takes the form: 

J= |^'d2/(l-2/){l^iv(v^2/) + ^^iv(v^[l-2/])}- (19) 




FIG. 2: Graphic illustration of the simplest term C. 

The similar fusion takes place in the Trotter integrals at an arbitrary order of ^-expansion. 

The further simplification is possible only in doing the thermodynamic limit N oo. In 
this limit the function Vn{ii) defined in Eq. lfTTjl collapses to the 5- function: 

P;v(v^y)U-.oo^^%). (20) 

Then the integral ()19|) is easily calculated and we find finally the expression for C: 

C = (6iV|r|) 7^7v(l), (21) 

where 71n{X) is given by Eq. lfTTjl . 

Let us look at the calculation of C from a different viewpoint. Namely, we note that in 
this calculation only two species of levels i and j are involved. We show this situation 
graphically in FigQ. Shadowed boxes mark the energy levels with different shadowing 
("colors") for Ei and Ej. They are connected by the interaction line (IL) which is associated 
with the factor [NrhY . Because of the fusion of two ©Ar-functions, the interaction of two 
levels brings only one 5-function with the normalization factor l/(A^|r|). As the result, we 
get one 'free' parameter h which is decoupled from the combination A^|r| that involves the 
matrix size N: 

c ~ ^^^7^;v(l) = X {N\T\hf 7^^(l) . 

It is clear from the above analysis that if we continue expanding in higher powers of Vij = V*^ 
but keep the number of colors equal to 2, we will increase the power of the combination 
(6A^|r|) but not the power of the free parameter h. 

C. Selection rule 

Consider a generic term of order IP'^ in the perturbation series: 

C2fe(A,B)= CN{{ki}) b"" {N\r\bf; A+B = 2fc. (22) 




FIG. 3: Example of diagrams classified in accordance with the selection rule: a) the case of two 
colors: k = ki = 2; h) and c) two simplest diagrams for the case of three colors: (b) k = 2,ki = 
1, A;2 = 1 and (c): k = 3, ki = 1, k2 = 1, ks = 1. 

It corresponds to a diagram with k = ki + k2 + ■ ■ ■ interaction lines distributed in a certain 
way {ki} = {fci,/c2,...} among the links connecting a given number of different energy 
levels {^^ colors"). 

Let us formulate the selection rule for such a perturbation theory. 

• The power A of the "free" parameter b is equal to the number of interacting energy 
levels minus 1: A = c — 1. 

• The sum A + B is equal to 2k . 

We exemplify it in FigOl for three different diagrams: the diagram with two colors (a) 
corresponding toA; = 2,A=l, B = 3; and two diagrams with three colors: (b) for 
k = 2, A = 2 , B = 2; and (c) for k = 3, A = 2 , B = 4. 

Suppose that we want to derive the perturbative term Kf^^r) of the order of 0{b^), 
which would be valid at an arbitrary value of A^|r|6 in the limit oo. Such physical 

result can be obtained using the following strategy: 

1 . Fix the power A of the "free" small parameter b by fixing the number of the interacting 
energy levels, i.e., the number of colors: c; 

2. Perform summation of an infinite series in A^|r|6 at fixed c over the number(s) ki of 
IL along each link. 

3. Analyze the limit ^ oo of the corresponding infinite series in A^|r|6. 



If this virial expansion works the series should converge in the hmit N oo. 

The strategy we suggested is a certain way of summation of perturbative series. Its 
physical meaning is analogous to the idea of resonant energy levels^^ or locator expansion^i. 
However, in contrast to these essentially heuristic approaches our selection rule offers a 
formalism that allows to implement a regular expansion in the small parameter b. It takes 
into account resonant as well as non-resonant levels accurately which is impossible in the 
above heuristic methods. 

Concluding this section we would like to note that the coefficient C^^lki}) in Eq.ip^ is 
a product of two different factors: 

CNi{h})=nNi{h})vi{h}), (23) 

where V{{ki}) is a universal coefficient which depends on a combinatorial factor /C and on 
the corresponding Trotter integral X. The factor 7l{{ki}) is not universal. It arises because 
of the summation of the product of the correlation functions — j\) over the real space 
and is a generalization of the function TZ^ik) given by Eq.()17|). In the large- limit it can 
be represented as: 

nN{{h}) = n{{h}) + l7^«({A;,,}) + . . . (24) 

where 7^({fci}) and VS^\{ki}) are of the same order. In this paper, we study mainly a 
contribution of the leading term TZ{{ki]). The role of the corrections to TZ{{ki}) is 

discussed below in Section fVII Bl 

IV. THE CASE OF TWO COLORS 

In this section we calculate the term of order h in K{r). According to a general scheme 
outlined in the previous section this correction is governed by the two-color diagrams: 
which correspond to the following power series: 

oo 

bki = Y.i-'^f ^(f^) T^Nik) ; Vik) = }C{k) I{k) ; x = iV|r|6. (25) 
fc=i 

The main difficulty in this calculation is the universal factor V{k) in Eq. ()25|l . Applying 
the Trotter formula to both traces in Eq.® and making an expansion in V analogous to 
Eq.()14j) we observe that the same number of IL may correspond to completely different 
distribution of the off-diagonal matrices V over two traces (see Fig.5). 
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FIG. 4: The two-color-diagrams. 

The combinatorial part of the calculation consists of two tasks. One of them is to dis- 
tribute colors (denoting Sij) over 2k segments. Since V is the off-diagonal matrix, the 
adjacent colors in each circle must be different. In general this is a particular case of the 
famous combinatorial problem of graph coloring^S. For the case of only two colors {si and 
Ej) considered in this section the problem is trivial: each circle must contain an even number 
of the segments while the first segment of the Trotter chain corresponds either to or to Ej 
(see Fig. 6), so that only the factor of 4 arises from this coloring. The second combinatorial 
task is to find the number of ways to connect V^-matrices in pairs making either in-circle or 
inter-circles connections. Such connections are equivalent to the Gaussian averaging. 

In Fig. 7 we show three different ways to connect two \^-matrices. Two links labeled by 1 
and 2 are in-circle connections, while the third one (labeled by 3) is the inter-circles connec- 
tion. The difference between the connections 1 and 2 is that the connection 1 corresponds 
to the average (VijVji) while the connection 2 corresponds to the average (VijVij). For real 
Hermitian matrices (the orthogonal ensemble, GOE, j3 = 1) Vij = Vjt both connections 




FIG. 5: Two examples of different distribution of y at k = 4. The first segment of the Trotter 
chain is denoted by the shadowed box. 



are possible and the corresponding averages are equal to each other. However, for complex 
Hermitian matrices with the same variance of the real and imaginary parts (the unitary 
ensemble, GUE, (3 = 2) only one of the two averages, {VijVji)^ is non-zero. In this case the 
"crossed" connection 2 is not allowed. 

Now one can easily find the combinatorial factors associated with the number of possible 
connections. In the GOE any of the 2k vertices V can be connected with any other V^-vertex 
and we get: 

1Cgoe{^)=^ x(2A;-1)!!. (26) 

In the GUE the 2k K-vertices are divided by two groups: one group containing k vertices of 
Vij and another group containing k vertices Vji. The connection is possible only if y- vertices 
belong to different groups. Thus we obtain for the corresponding combinatorial factor: 



ICcuEik) = 4 X fc! . 



(27) 



In order to find the coefficient T{k) we have to fix the number of V matrices in each 
trace to be 2m and 2{k — m), thereby also fixing the number of the segments in the circles. 
Each segment corresponds to a Trotter variable to be integrated over. The details of the 
calculation of the Trotter integral T{k) are presented in the Appendix |B| The final result 
is obtained after summation over m: 



I{k) = b - 



1 



1 



N\T\b 2 k\{k- 1)1' 
Substituting the combinatorial factor /C(/c), Eqs. ()26H?r|) . and the Trotter integral 



(28) 





FIG. 6: Two different ways to put colors in the two-color problem. The first segment is denoted 
by the shadowed box. 



FIG. 7: An example for 3 different connections of ^-vertices at A; = 4. 
into the series we find: 

oo 

k, = l^f} Y.i-^fcf{k)nr,{k) {N\r\hf'^-'- (29) 



fe=i 



= jj^, ■ (31) 



Eqs. fl29llt?T|l correspond to the leading in 6 ^ 1 correction to the Poisson spectral form- 
factor. They have been obtained in the framework of the two color approximation when 
we have taken into account only pairs of the interacting energy levels. We note that the 
transformation (j20|l of the DAT-function into the delta-function is approximate and its sub- 
leading terms yield the corrections to hKi of higher orders starting from ~ })^{^Ki). 

We emphasize the the formula for Ki is valid for a Gaussian random matrix ensemble 
with a generic variance of the off-diagonal matrix elements = h"^ — 



V. THE CASE OF THREE COLORS 

According to the selection rule explained above, the term of order 6^ in K{t) results from 
an interaction of three energy levels, i.e., we have to perform a summation of the three-color 
diagrams, see FiglHl Three shadowed boxes stand for the three independent energy levels 
Ei, Ej, and El {I > j > i), which are connected by ILs. Any two levels can be connected 
with each other independently of a connection with the third one by an arbitrary number 
of ILs: ki, /c2, and k^. This leads to a three-dimensional power series in a; = A^|r|6: 

oo 

b'k2 = E Viih}) TlNiih}) ; 2 = 1,2,3. (32) 

^1,2,3 



Here, the universal factor V{{ki}) contains the Trotter integral 2{{ki}) and the combina- 
torial factor IC{{ki}). The function T^jv is a generalization of the sum in real space for the 
case of three colors: 



N N N 



nN{{h}) = ^ E E E 

l=j+l j=i+l 1=1 



(33) 



Generically, any energy level interacts with the two others and all numbers of ILs are not 
zero fci,2,3 > 0. A diagram describing this case will be referred to as a triangle (see FigJHlA). 
However, there are configurations where two levels do not interact with each other and one 
parameter ki is zero. A corresponding diagram will be named a tine (see FiglHtB). We can 
attribute a certain physical meaning to both of the three-color diagrams. Let us involve a 
tree-like structure in the real space where a path with loops is prohibited. It will generate 
only the "lines". On the contrary, a real space structure where the paths with loops are 
allowed will yield both the "triangles" and the "lines" . 

The combinatorial factor JC includes the number of coloring of segments by three colors 
(denoting Sijj) and the number of ways to pair ^-matrices by in-circle and inter-circles 
connections. The V^-matrices create three sets of vertices in the circles: 2ki vertices of 
{Vi,i,yi,i}y 2A;2 vertices of {V^-,;, VJj}, and 2k3 vertices of {Vij,Vj^i}. In view of the property 

{VsM = (V^i) Ss,, 5p,r + {\Vs,,\') 6s,r (34) 

the vertices belonging to one and the same set must be paired with each other independently 
of the two other sets. Applying the same arguments that were used in the case of two colors 




FIG. 8: The three-color-diagrams: a) a triangle diagram with ^1,2,3 > 0; b) a line diagram with 
^2,3 > and ki = 0. 



we arrive at: 



ICgoe = CoIgoe n (2fc.-l)!!; (35) 

j=l,2,3 

K^GUE = CoIgue n (^i)'? (36) 

j=l,2,3 

where Col is the coloring factor. When the number of colors is larger than two the coloring 
factor is different in GOE and GUE. Details of calculation of Col are given in Appendix O 
Here we explain notations and give the results for GOE and GUE. 

Denote a total number of the segments for each of the three colors by R, G, and B. 
These numbers can be expressed in terms of ki^2,3- 

R = kl + k2; G = k2 + k3; B = kS + kl . 

The colored segments are distributed over the two circles. The numbers of the colored 
segments in the first circle will be labeled r, g, and b. Thus the second circle contains 
R — r, G — g, and B — b segments of different colors. 

In the case of GOE, the two circles are colored independently under the only condition 
that any two neighboring segments must be of the different colors (see Appendix E))- The 
factor CoIgoe reads 

CoIgoe = CoP^ (r, g, b) x CoP^ {R - r,G - g, B - b) ; (37) 

a-\-b~\-c 

Col'^^^ (a, b,c) = {a + b + c) 2^+^+^ V — ^ ,,7' , — . 

.=n.atra,M - «)' (/^ - &)! (/^ - c)!(a + & + c - 2/.)! 

Before exploring the case of GUE, we remind that the "crossed" connections (see example 
2 in FigEj) are not allowed in GUE. The "crossed" connections can be avoided if each set 
of the vertices {Vp^q, Vq^p} contains exactly ki of Vp^g vertices and ki of Vq^p vertices 
regardless to their distribution over the circles. Note, that this restriction applies only to 
the total number of Vp^q and Vq^p while the number of the vertices Vp^q and Vq^p in the 
same circle is not necessarily balanced. The coloring of the two circles is not independent 
any longer: the difference d = Number [V^^^] — Number [V^^p] of the conjugated vertices must 
be the same in magnitude (but opposite in sign) for the two circles. Therefore one has to 
compute the coloring factors for each ring at a fixed difference d and then to perform a 
summation of the product of two coloring factors over d (or, equivalently, over the related 



variable D = Mr + b + g— \d\)). The result reads (see Appendix IHl for the details): 



ColcuEir, 9, b, R, G, B) = hk E ^D) x 

D=max(r,(/,6) 

I j2 {h - D - P, D - P, P - r, p - g, p - b) ^ ^^^^ 

vp=max(r,g,fe) -P 



E 



(^i2M^_D-q/-^ + D-q,q-{R-r),q-iG-g),q-{B-b) 



g=msLx{Fi,~r,G—g,B—b) ^ J 

h^r + g + b, l2 = R-r + G-g + B-b; 
\ 2, ifZ}<^; 

1 , if D = — and — is integer . 
(ai + 0:2 + tts + . . .)! 

Brackets (ai, 0^2, 0^2, . . .) = ; — ; — ; denote the multinomial coefficient. 

ai! 02' as' ■ ■ ■ 

The coloring factor has much simpler form for the line-diagram. For the example shown 
in FiglHlb, where g = r + b; G = R + B, Gol acquires the following form in GOE and GUE: 

Gol\g=r+b-^G=R+B = Coliine{r, b) Golline{R - r, B - b) ] (39) 

{a + b)\ 



Coliine{a, b) = 2- 



a! 6! 



There is no essential difference in the calculation of the Trotter integral for the case of 
two and three colors. We follow the algorithm explained in Appendix El For three colors, 
the integrand contains a product of two "infusible" T>n functions. They are converted to 
the 5-functions and the result of integration with the accuracy of l/(A^r)^ is 

' ' ' {Nry hh {R + G + B-4y. ^ ^ 

1 1 1 

X ■ 



(r -l)\{R-r- 1)! {g - 1)\ {G - g - 1)! (6 - 1)\ {B - b - 1)! ' 

where the sum over all possible arrangements of the first Trotter segment is done. 

To derive the universal factor V{ki, k2, k^) in Eq. ()32|) . we have to perform summation 
over r,g, and b accounting for the different distributions of K-matrices over the traces. We 
remind that the matrices V from the ffist trace enter the Trotter chain with the minus 
sign. Therefore, the summand must be multiplied by the factor (— 1)'^. Unlike the case of 
two colors, both the coloring factor ()37| IHH|l and the Trotter integral pn|l depend on the 



distribution of the colored segments. Thus we arrive at the sum: 

R~-\ G—\ B 1 

r=l 3=1 b=l 

which turns out to be rather cumbersome. Amazingly, these summations can be done 
analytically and we finally find the correction K2 to the spectral form-factor: 



^ fci,fc2,fc3=0 



(41) 



where TlN^ki, ^2, k^) is given by Eq.(33), x = A^|r|6 and the coefficients are 



^(3) ^ r(/ci + fca + ks) H^=i(fci) H^=i(fc2) S^=i(fe3j ..^s 

T{k, + k2 + k3-3/2)T{k, + k2)T{k2 + k3)T{k, + k3)' ^ ' 

2'=r(A;-l/2)r(A; + l/2) 



S;3=l(fc) 



^V{k + 1) 



^/3=2 



r(A:- 1/2) 



(fc) = 



r(fc + i) 



The correction is governed by the interaction of the three energy levels and results 
from the summation of "triangles" and "lines" : All terms with ^^1,2,3 > 1 correspond to the 
"triangle" diagrams while the "line" diagrams are obtained by setting /ci = or k2 = or 
A;3 = 0. 

The series in r.h.s. of Eq. Mlj) is three-dimensional and cannot be reduced to a product 
of one-dimensional series. Therefore an analysis of its behavior at arbitrary x is not trivial. 
In the case of GUE, one can represent the function r~^(A;i + k2 + k^ — 3/2) as an integral 
using the identitySi; 

1 1 /■°° exp(a + it) , „ 

' , a > , (44) 



T{z) 27ri J-00 (a + ity 
and change the order of summations over ki and integration over t. The real-space summa- 
tion which is implied in the function TZn has to be done at the last step. This is an effective 
tool for the asymptotic analysis of the series ()4ip43|l . The case of GOE is more complicated. 
Numerical and semi-numerical methods are useful here. We will demonstrate an application 
of Eqs. ()4HI43|) to the Rosenzweig-Porter model and to a critical almost diagonal PLBRM in 
the forthcoming paper—. 



VI. THE CASE OF CROSSOVER: ALMOST UNITARY GAUSSIAN ENSEMBLE 



We end the presentation of the method with the consideration of an almost unitary 
Gaussian RM ensemble: we explore the case of a crossover between GOE and GUEj^ which 
is close to GUE. To define the almost unitary ensemble we introduce a parameter rj into 
Eq.()34j). which controls an asymmetry for the variance of real and imaginary parts of the 
hopping elements: 

{Vs,pVg^r) = b'^J^i\s-p\) [ 5s,rSp,g + r]Ss,q5p,r ] ■ (45) 

In GOE, Vs^p is a real number and tjgoe = 1- hi GUE the entries Vg^p are complex 
numbers, where 3?[K,p] and are statistically independent random variables with 

equal variances and therefore rjcuE = 0. Hence the parameter f3 is implicitly linked with 
rj. An ensemble of RM is almost unitary if < r/ <^ 1. 

A first (trivial) impact of small t] on the form-factor arises from the proportionality: 
Ki ~ ^/J3 ; K2 ^ P (see Eqs. fl^ 1^ ). which comes from the Trotter integrals. We will 
concentrate on the other outcome of a deviation from GUE where the non-zero parameter 
rj leads to an appearance of new series in powers of x. We will derive a series that yields a 
leading in r] correction to -ft'i,2|/3=2- 

When one performs the Gaussian averaging of ^/-matrices, l^-vertices in the circles can 
be paired by either "crossed" or "non-crossed" connections (see FigEj). Each "crossed" 
connection is proportional to rj in accordance to the definition ()45p . Therefore, the "crossed" 
connections are prohibited in GUE ( r/ = ). To explore the leading term in 77 -C 1 in the 
almost unitary case we have to select the diagrams with a minimal number Cr of "crossed" 
ILs. Our main task is to calculate their combinatorics. 

We start with the two colors. Let us consider a given distribution of ^/-matrices over the 
traces at the fixed parameter k in GUE. There are sets of conjugated vertices: k vertices 
Vij and k vertices Vjf, and there are k connections without crossing. Each "non-crossed" 
connection pairs the vertex Vij with its conjugated counterpart Vji . A diagram with the 
minimal number of "crossed" ILs is obtained if we choose two "non-crossed" connections 
and rearrange them into two "crossed" connections, see an example in FiglHl The rest 
^ - 2 ILs remain "non-crossed". There are {kik - l)/2)^ ways to choose two V.,- and 
two Vj-^j- vertices for the two "crossed" connections and (A; — 2)! ways to draw the remaining 



FIG. 9: Substitution of the "non-crossed" connections by the "crossed" ones: In GUE, only the 
dotted "non-crossed" hues are allowed; the diagram with the minimal number of the "crossed" 
lines is obtained if one substitutes two dotted ILs by the solid ones. 

"non-crossed" ILs. Therefore the combinatorial factor for the diagram with the minimal 
number Cr = 2 reads: 

;C,«i(A:)=4x^^i^^ A:!, (46) 

where the factor 4 accounts for two different ways to color each circle by two colors. Sub- 
stituting Eq. ()46j) and the Trotter integral p8|l with [3 = 2 into the two-color series (j25|l . we 
find the correction: 

AKi = ^V' E(-l)' (^3^ T^nW x"^-' ■ x^N\T\b. (47) 

Each three-colored diagram in GUE has 3 sets of conjugated vertices, which numbers 
are balanced: Number [V^^;] = Number [VJ^ J = ki, Number [V^-;] = Number [V^j] = ^2, and 
Number [^j] = Number [V^-^j] = k^. To generate the diagram with the minimal Cr = 2 we 
can repeat the arguments explained above for the case of the two colors: a pair of "non- 
crossed" connections can be rearranged into a pair of "crossed" ones successively in the 
bonds Vi^i ^ Vi^i ; Vj^i ^ Vij ; and Vij ^ Vj^i. This procedure results in the combinatorial 
factor 

^ 1=1,2,3 

with the coloring factor CoIgue defined by Eq. ()38|l . Therefore the first correction to GUE 
takes the form: 

fcl,fc2,fc3=0 



^(3) _ h ih -l) + k2 {k2 -l) + h ih - 1) (3) 

Another way to insert the "crossed" ILs into the three-colored diagram in GUE is to 
violate the balance of the vertex numbers Number [Kj^;] — Number [VJ^j] = di, Number [V^y] — 
Number[V;j] = d2, Number[Vij] — Number[Vj- j] = d^, di^2.3 = ilj =t2 . . ., and to compensate 
this violation by the "crossed" connections. One can prove that the integer parameters d^ 
are always equal to each other: di = d2 = d^ = d (see Appendix IH}. Therefore the minimal 
number of the "crossed" connections obtained in this fashion is Cr\d=±i = 3 which is, 
however, larger than Cr\d=o = 2 in Eqs. ()49ll!Tn|l . 



VII. EXAMPLES 



In the last section, we give two examples of the application of the general Eqs. ()29ll^ 
im imilSUj) for the corrections bKi and b'^K2 to the Poissonian form-factor Kp = 1. 
The examples are the Rosenzweig-Porter model with the parameter b depending on the 
matrix size and the power law banded RM with the size-independent parameter b. 



A. The Rosenzweig Porter model 

Let us consider the Rosenzweig-Porter model with 

6=|,i3«l; T{z-j) = l. (51) 

This definition corresponds to an almost diagonal RPRM where a weak perturbation of the 
diagonal matrix can nevertheless yield a nontrivial level statistics^^. 

We note first, that the parameter x does not depend on the matrix size A^ : 

N 

x = N\T\b= —\t\B^T. 
II ^11 

The ratio N/N is a constant of the order of 1. The product b'^^^TlN{{ki}) ( c is the number 
of colors) also has the finite limit at A^ ^ oo. For example, in the case of two colors 

hm {bUNik)] = lim ( — ) = - ; 

and Eq.(j2ni) takes the form: 

no 

\k /^{2)/7 \ rp2k-l 



bK,U^^ = ^nP B Y^(-lfc';\k) T''-\ (52) 
fc=i 



(2,) 

with the coefficients {k) from Eqs 

Let us consider for simphcity the complex hopping elements, i.e., (3 = 2. The sum over 
k in Eq.(jS21) can be easily calculated and we find: 



hK^ {T,(3 = 2)\N^^ = -V2n BTe-^ . (53) 



In a similar manner we treat the case of three colors, where 

lim (6^7^A^(A;l, ^2, k^)) = — ; 



and consequently 

b'K^lN^oo = ^S' ii-^r C,is) T'^-' . (54) 



=2 

The coefficient C(^{s) is given by the sum: 



3 J 



kl,k2,k3=0 s'=l fc3=0 



Here we have introduced the sums of two and three indices: s = ki + k2 + k^ , s' = k2 + k^; 
and is taken from Eqs. ()42|4Hp . In the case of the complex hopping elements, /5 = 2, 

all summations over k^, s', and s in Eq.()54|) can be done analytically. The correction to the 
form-factor takes a simple form: 

b'K2iT,(3 = 2)U^,o = i2T' - 3) e"^^ . (55) 

y 

The results fj53p and (jSHl) are in agreement with the expansion in powers of B of the exact 
expression for the form factor K{T) obtained for /3 = 2 in the paper [lo| . 

The case of the real hopping elements, (3 = 1, can also be analyzed based on the formulae 
but the calculations are lengthy and we present them in the forthcoming paper. 



B. Critical PLBRM 



Before giving the example of a critical PLBRM, we would like to discuss generic properties 
of the virial expansion in the case of the iV-independent parameter h and a converging real 
space sum Eq. (II 7|1 . In this case, all corrections to the form factor Km are given by the 
series in powers of the parameter x = A^|r|6 that diverges in the thermodynamic limit. 
In the leading in approximation in Eq. p^ . where 'R-N{{ki\) = 'R-{{ki\), the entire 



dependence of the spectral form-factor K{t) on r comes only through the dependence on 
the parameter x. This means that K{t) is either r-independent or divergent (if the series 
in A^|r|6 determines a function that does not have a finite limit as A^|r|6 ^ oo ) in the 
thermodynamic limit. 

In order to obtain a dependence on r for a finite limiting spectral form factor /r(r)|jv-^oo 
one has to account for 1/A^-corrections to the real space sum 7l{{ki}) in Eq.()24|). Rewriting 
~ b\T\/ (A^|r|6) and absorbing the factor (iV|r|6) from the denominator into the infinite 
series, we obtain the correction of order br to the limiting spectral form factor: 

^(r)liv^oo ~ Xo + Xib\r\ + ... 

where Xo(^) = 1 + Cqi 6 + C02 + • • • is the unfolded level compressibility (see Eq.((7j)); 
Xi{b) = Cub + cub"^ + . . ., and Cij are numerical coefficients. Thus the correction to 

7l{{ki}) determines xi ^ind leads to the 6^|r| correction in the limiting form factor. For 
6^1 this correction is small even at the Heisenberg time r ~ 1. However, it could be 
important for the tail of the two- level correlation function Eq. Q at large energy separations 
uj because it is non-analytic in r. 

Now, we illustrate an application of the method developed using an example of a critical 
almost diagonal PLBRMi^: 

6 = S«1; J^(\^-J\)^l-^—^, ^^J. (56) 

We restrict ourselves to the case Br <^ 1 omitting 1/A^-corrections to the real-space 
sum. The leading term of TZN{k) takes the following form: 

1 ^ 1 

n{k) = -r lim y (57) 

m=l 

Let us substitute expression ^7\i into Eq. (j^ and explore an asjnnptotic behavior x = 
N\t\B ^00 of the series pUHnj) and (MTTIl for GOE and GUE respectively. The simplest 
route is to change the order of the summations: 

The summation over k yields: 

. £ 4 -P f-^) I - • " ^ ^ ^ (59) 

N\t\ ^^w? \ Im^j 1^ ^ = 2. 



Here /o,i(' ■ ■) are the Bessel functions. The sum over m converges at m ~ x ^ 1 therefore 
it can be converted to the integral over m. After this integration we find the two-colors 
correction to the level compressibility xo — 1 + x coi: 

Coi 1/3=1 = -2 ; Coi 1/3=2 = -TT . (60) 

Note, that Eq.^ can be reproduced with the help of a heuristic renormalization-group 
method (compare with [16]). 

Using our method it is also easy to derive from Eq. ljTfj) the correction in the almost 
unitary case: 

We emphasize that it is the behavior of J-'{\i — j|) at large |^ — j| that affects the factor 
TZ{k) and governs the asymptotic behavior of the corresponding series. Depending on it, the 
two-color correction may either approach a constant limit or diverge as ^ oo. One can 
show^that: 1) for J-'{\i—j\) decreasing /aster than l/|z— the two-color correction vanishes 
in the limit — * oo; 2) for jF(|i — j|) decreasing slower than — it is divergent. This 
is the manifestation of the localization-delocalization transition in the spectral statistics. 
Note that even the case, where the series in r.h.s. of Eq.()29|) does not have a finite limit as 
A^ oo, can be considered by our formalism as long as the correction Ki{t) is small at a 
finite A^. 

A remarkable feature of the critical PLBRM with JF ~ l/\i — in the three-color 
approximation is the logarithmic divergence log^(A^|r|i3) of the line- and the triangle- 
diagrams (see FiglH)) at any (3. However, in the sum of the line- and the triangle-diagrams, 
the logarithmically divergent terms of order cancel out for both GOE and GUE so 
that \imx^oo(^K2{(3 = 1, 2)^ is a finite constant of order 1. The sub-leading diverging term 
survives strikingly in the three colors correction of the almost unitary ensemble: B^AK2 ~ 
—B'^rj'^ log{N\T\B). This indicates a failure of the virial expansion for the critical PLBRM 
as A^ — s> oo in the crossover between the unitary and orthogonal ensembles and raises a 
question on the cancellation of the logarithmically divergent corrections of higher orders in 
& in a pure GOE or GUE. This question deserves a separate detailed study elsewhere. 



VIII. CONCLUSIONS 



In this paper, we have developed a new method to study the spectral statistics of Her- 
mitian Gaussian random matrices with parametrically small hopping elements ifjj as com- 
pared to the diagonal ones: {\Hi^j\'^) / {e\) IP' <^ \. We have derived a regular virial 
expansion of the spectral form-factor in powers of 6, where the virial coefficient in front of 
h-^ is given by interaction of ^ + 1 levels. 

The expansion is represented by diagrams which are generated with the help of the Trotter 
formula. We have established the rigorous selection rule for the diagrams, which allows to 
account for exact contributions of a given number of resonant and non-resonant interacting 
levels. Thus the method offers a controllable way to find an answer to the question when 
a weak interaction of levels can drive the system from localization toward criticality and 
delocalization. 

The method applies to the spectral properties of the random matrices with uncorrelated 
entries and with a generic dependence of the variance = 6^jF(|z— j|) on the distance 

\i — j\ from the main diagonal. We calculated the corrections governed by interaction of 
two, Eqs. ()29li;TT|l . and three, Eqs. ()4Hl4H|) . levels in the cases of GOE, GUE and in the almost 
unitary ensemble, Eqs. ()47l49ll50p . for an arbitrary function J^{\i— These equations are 
the main results of the paper. We have demonstrated their application to an almost diagonal 
Rosenzweig-Porter model and to an almost diagonal power law banded random matrices. 

The calculation of the virial coefficients is based on the solution to the combinatorial 
problem of the graphs coloring. The solution for the coloring of a single graph is known 
in a closed form for an arbitrary number of colors^^. Thus using our method one can, in 
principle, compute the virial coefficient of an arbitrary order. 
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APPENDIX A: THE LEVEL COMPRESSIBILITY 



In this Appendix, we derive an approximate relation between the level compressibility 
and the correlation function R. Let us define the two-level correlation function Re{}jj) 
centered at an energy E: 

o r ^ {{i:n,m5{E-e^)5{E-^-e„.))) 

^^(^) = (Al) 

In this definition it is assumed that u is small compared to the total bandwidth, so that 
{p{E))^{p{E-uj)). 

We set apart the diagonal (singular) part in the double sum in numerator of ()A1|) : 

i?^M = i?J)M + ^|||y. (A2) 

Then the normalization by (p(-E))^ ensures that the regularpart R^p{uj) depends very weakly 
on the energy E, except for regions of the width ~ u near the band edges. 

The level number variance can be expressed in terms of the correlation function R{u) j^]: 

^2{n) = J^ i, d^dE'. (A3) 

Substituting decomposition ()A2jl into Eq. ()A3|l . one finds 

S,(n) ^ n + J2 i?«(.) ^ ; i?(^)(.) ^ 41o(-) • (A4) 

This expression is valid up to the leading order 0{b) with respect to 6^1. 

The definition (jSj) for the level compressibility yields now the following expression: 

Y = 1 + lim r ds lim R'^^'^sA) . (A5) 

Using definition (Al), formula ()A2jl and the identity 

/+00 r+oo 
dE / dcu((p(E)p(E-cu)))e'"* = ((Tre'^* Tre-'^*)) 
-OO J —OCl 



/+00 
dse'^'^R'-'-^sA). (A6) 
-OO 



-OO J — OO 

we derive a useful relation: 

1 , . /- + 00 

iv' 

which remains valid also in the limit ^ oo. 

At t = to = l/{nA), large energy scale u = sA ^ nA does not essentially contribute 
to the integral in Eq. ()A6|) due to the strong oscillations of the integrand. Hence, we have: 

lim r ds lim R'^^'^sA) = lim /^°°dse'"^° lim R^''\sA); tq = toA , 

and using Eqs.(A5,A6) arrive at the expression (|7|). 



APPENDIX B: INTEGRATION OVER THE TROTTER VARIABLES 



In general, the Trotter integral depends on the number of colors and on the distribution 
of matrices V over the traces. The former defines the number of 'infusible' functions 
in the integrand while the number of the Trotter variables is determined by the latter. In 
this Appendix, we calculate the integral in the case of 2 colors at an arbitrary number of the 
colored segments in the circles determined by the parameter m. The integrand in Eq. ()19|) 
(see the Section IIIIB 2|) contains two additive parts reflecting two different possibilities 
to color the first segments. Unlike that example, we will use another procedure with a 
simplification that we will not pay an attention to a specific color of the first segments. The 
final result of the calculation is not affected by this simplification. 

Let us apply the Trotter formula to both traces in r.h.s. of Eq.® involving two infinitely 
long Trotter chains of the length ni 2 00 and let the number of the colored segments 
in the circles be 2m and 2{k — m). We remind: correct coloring by 2 colors imposes a 
restriction that a circle must contain equal number of sectors for each color. Therefore, the 
total number of the segments for a given color is k with m segments belonging to the first 
circle and k — m segments in the second one. 

Denote the discrete Trotter variables coming from the first trace as yi, 1/2, ■ ■ ■ , y2m so 
that the variables yi, y2, ■ ■ ■ , Vm mark the segments of the first color while the rest of 
variables are attributed to the segments of the second color. The variables from the second 
trace are denoted as zi, Z2, . . . Zk-m for the first color and Zk-m+i, Zk-m+2, ■ ■ ■ Z2{k-m) for 
the second one. Since two neighboring subsegments separated by the origin of the Trotter 
chain (instead of the vertex) are of the same color we fuse them and mark by the single 
Trotter variable. 

We rescale the variable introducing the continuous ones: 1^ = yi/rii , Zi = Zi/n2- They 
are normalized as (see Eq. ()15|) ) 
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The integrand must be multiplied by the factor: 

2m 2{k-m) ^ ' 

The numerator reflects an invariance with respect to a cyclic permutation of the matrices un- 
der the traces while the denominator excludes multiple counting of the equal configurations 



of V-vertices in the circle. 

We follow further the algorithm explained in the Section IIIIB 21 and, accounting for the 
property ()B1|1 and fusing two Vn functions as it is done in Eq. lfTIUl . the Trotter integral 
can be presented as: 
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Having integrated over all variables but Y = Yi', Z = Yl'jZT '^h arrive at 
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At the last step, we sum over the parameter m to account for different distributions of 

\^-matrices over the traces: 

I(,).|j„.,„).v|£j__L_. (B6) 

APPENDIX C: COMBINATORIAL COEFFICIENT OF COLORING 

To calculate the universal factor V{{ki\) at any number of colors c we need a combi- 
natorial factor of the circles coloring: ^-vertices divide the circles into the segments. The 
colors, which denote the statistically independent energy levels e^, should be distributed 
over the segments by all possible ways under the condition that all adjacent segments must 
bear different colors. 

Let us start with the coloring of a single circle. The question of our interest was formulated 
in Appendix of paper js^ as follows: How many permutations of L objects of c different 
colors (given gi objects of the first color, g2 objects of the second color, etc.) are there 
under the condition that no objects of the same color may stand next to each other? We 
have to supplement it with two additions: 1) we do not distinguish the segments of the same 
color; 2) we can do cyclic permutations of the segments in the circle. Thereafter, the generic 
coloring factor can be obtained easily from Eq.(30) of paper j3]: 
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FIG. 10: a) An unrolled circle with r = 2 and with two clusters of green and blue segments; b) A 
cluster of odd length; c) A clusters of even length. Pairs of bold letters denote switching of colors 
at the y- vertices. 
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Eqs. ()CHIC2|) describe any situation with the arbitrary number of colors and the arbitrary 
number of the segments. Unfortunately, these expressions are cumbersome and are not 
convenient for a further calculation of the factor V. In this Appendix, we will derive a 
simpler formula for the particular case of three colors, c = 3. 

We use the notations from SectiorfVl labeling the numbers of the colored segments by 
f = 9ij 9 = 92, and b = g^. The first, the second and the third colors will be denoted as 
"red", "green" and "blue". Let a color of the first segment be, for example, red. There are r 
red segments which separate out r clusters. The clusters contain green and blue segments, 
which must be arranged in an alternating fashion. A number of the segments in the cluster 
will be referred as a "length of the cluster". The are different types of the clusters (see 
Fig lTn|l : 1) The clusters with an even length that enclose pairs of sectors; 2) The clusters 
with an odd length that enclose the pairs as well as one extra segment (either green or blue). 
We denote a number of these clusters accordingly: uq, Ug, and n^. The parameters no,c,,6 
are connected to the numbers of the segments: 

Hg + riQ + Tib = r ; g - b = Ug - rib . (C3) 

First we have to distribute the clusters of three types over r positions. The number of 
the distributions is given by the multinomial coefficient: 

(no, Ug, Ub) = —r— — 7 . (C4) 
nolnglnbl 



Given that each of no clusters of even length contains at least one pair we distribute 
b — rib — uq = g — Ug — riQ pairs of the segments over r clusters. The corresponding combi- 
natorial factor is: 



where 2"° reflects the possibility to start the clusters of even length with either a green 
segment or with a blue one. 

The coloring factor Colr{r, g, b) for one circle with the red first segment is obtained after 
summation of the product of Eq. ()C4|) and Eq. ()C5|) over the independent number of the 



second Colg{r, g, b) and the third Colb{r, g, b) parts, we have to repeat the same procedure 
assigning green and blue colors to the first segment of the circle. Eq. fl37|l has been obtained 
after the summation of all parta^: 



Eq. ljHHj) . that is used to calculate the universal factor V{ki, ^2, k^) at c = 3 in GUE, has 
been obtained under an additional condition. Namely, to exclude the "crossed" connections 
of the l^-vertices each set must contain an equal number of the "conjugated" 

vertices: Number = Number [V^^p] (see Sectional). This restriction applies only to the 
total number of Vp^q and Vq^p while the number of the conjugated vertices in the same circle 
is not necessarily balanced. Therefore, coloring of the two circles is not independent. 

Let us consider the balance between ^ and Vq^p in more detail. The conjugated vertices 
inside a cluster with an odd length and on its boundary are always balanced (see an example 
in FiglTUb). Hence, the clusters of odd length cannot violate the total balance. It is not true 
for a cluster of an even length: the imbalance between the conjugated vertices is ±1 for all 
three sets: {V^,z,V;,i}; {Vj^uVij}; and {Vij.Vj^j} (see an example in FiglTnh). The sigh of 
the imbalance depends on a color of the first segment of the given cluster. 

Denote the number of clusters of an even length, which start with green and blue seg- 
ments, by no~^ and Uq"'' = nQ—no'^ respectively and introduce a difference d = u'q'^ —nQ~^ = 
nQ — 2nQ~^ which will be referred to as "the number of defects". It gives the number of imbal- 
anced conjugated vertices of each set in a given circle. Below, the upper indices in and 
d^"^^ will label the number of defects in the first and in the second circles correspondingly. 




(C5) 



clusters, for example, over the parameter Ub in the range < n?, < 



r+b—g 
2 



To find the 



Co/(=^) (r, g, b) = Colrir, g, b) + Colg{r, g, b) + Colb{r, g, b) . 



Since each cluster of even length produces the same imbalance between Vp^g and Vg^p in 
all sets of the vertices, the number of the conjugated vertices is balanced if 



(C6) 



Equality ()C6|1 must hold true at any combination of the colors of the first segments and 
expresses the difference of the coloring in GOE and GUE. Eq. ()37|l . which is used in GOE, 
results from an independent summation over the number of defects di and d2 in each circle. 
In GUE, we calculate the coloring factor Coir of one circle at the fixed d and perform the 
summation over the number of defects obeying Eq. ()G6|) : 



J2 E { Colrir, g, 6, X ColriR -r,G-g,B-b, S^^) x } . 

The summand is the product of the coloring factors for the two circles. 

To find the expression for Coir we have to sum the product of Eqs. ()G4|) and the factor 
similar to Eq. ()G5|) : 
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over the parameter n^. The factor 



n 



^0 1 . np-d 

with ' = = 0, 1, 
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accounts for the 



different distributions of the clusters of the even length starting with green and blue colors. 
Eq. ()H8j) results from the summation over all combinations of the colors of the first segments. 
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